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[57] ABSTRACT 

A method for determining the coverage of a pixel, 
which includes determining a function of the distance 
from the currently sampled point to each edge of a 
polygon and then adding a predetermined value to the 
value of the function. This value is then clamped ac- 
cording to a function which provides an output of a 
predetermined value if the input to the function is not 
within a predetermined range and otherwise provides 
the value of the input if the input is within the predeter- 
mined range. The coverage of the pixel by the polygon 
is computed by multiplying the clamped functions. 
Blending the parameter values is performed if the value 
of the coverage falls within a predetermined range; 
otherwise, no blending occurs. 
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„ and the integer X values are shown as columns. A pixel 

METHOD FOR DISPLAY RENDERING BY center is assumed to be at the intersection of each row 

DETERMINING THE COVERAGE OF PIXELS IN and column, for example as shown at X= 11 and Y=4. 

POLYGONS It will be appreciated that each column is parallel to an 

5 orthogonal axis of the display device. 
r»AK?n£ ci ^°^^L a PP Ucation Sen No " In the prior art methods for scan converting, parame- 

07/657,087, filed Feb. 19, 1991. ter values are computed at both edges of a triangular 

BACKGROUND OF THE INVENTION polygon which intersect a particular scan line from the 

1 u- u C4X. t parameter values at the appropriate vertices. For exam- 
L. Field of the Invention 10 p i e , for scan line Y=«, parameter values along the edee 

JSSS^S 1 ^ ^ ? ^ ™ V2 be at location 12 (wWch is nofa 

xng images in a computer controlled raster dispkvy sys- pixe l center) and the parameter values for location 15 

T^ZZl^^i the present mven0 °? rela ^ 10 along edge V1-V2 will also be computed tS vies 

an apparatus and method for scan converting shaded ^ ±e ^ of ^ , ^ * uted b m 

SSnSf.r rdet 7^^«Wofp«fc» 15 ktinl b^tweln the palter vduesatongtheSu- 

pdygons m a computer controUed raster display sys- ^ edge. For M u ^ 

2 Prior Art P 3 * 310 ^ values computed by interpolating the 

Manycomputersystemsto^^ ^2^1^ ^iX? ™ " ™f ^ 
ing complex 3-dimensional curves and solids on display 20 ^^tf f 8 ^ ^f*"* aU 
devices which are controlled by these computer sys- mtersectl0ns each edge of the polygon with aU scan 
terns. These complex 3-dimensional solids and curves . ruiming through the polygon, then actual parame- 
are typically displayed interactively to allow the user to ter ™ ucs ?° r oenters m ±cn interpolated along 
manipulate the object Well known graphics techniques , scan . by mter Polating between the parameter 
for rendering these 3-dmensional curves and solids, 25 valu « at Ae edges along that scan lme. Thus for exam- 
such as Gouraud shading, hidden surface algorithms, P 7 , P arain eter value for the pixel at 8,5 will be 
clipping, filling polygons and coordinate transforma- c^culated by mterpolating the parameter values at loca- 
tions are used in order to generate the displayed object Uo T n 11 ^ d 16 ? f FIG * 1 

on a display device, such as a CRT video display which 1x1 811 aJtemative embodiment of the prior art method 

is controlled by the computer system. 30 of sx ^ n converting filled polygons, parameter values at 

These various computer systems often rely on con- ^ intersections of each column with each edge are 

verting a polygon representation of the solid to a raster calculated and these values are used for interpolating 

scanned image which is stored in the frame buffer of the parameter values at pixel centers along each column, 

computer system. In this circumstance, it will typically Thus, parameter values at locations 21 and 22 would be 

be necessary to convert from the vertices of the poly- 35 calculated by interpolating from the parameter values at 

gon which are calculated in the clipping and transfer- ^ V2 for location 21 and from the parameter 

mation manipulations to an array of parameter data values at V0 and VI for location 22 and then the param- 

which is stored in a frame buffer array at the resolution eter va3ue at P^el center X=3, Y=4 would then be 

of the display device. In other words, it is necessary to calculated by interpolating between the parameter val- 

convert from the specification of parameter values (e.g. 40 ues at locations 21 and 22. 

color in the form of "RGB" or other representations of * l ™& De appreciated that considerable computa- 
color and Z values and other useful parameters) at each tkmal effort is required in order to first calculate param- 
vertex of a polygon to parameter values at each picture eter values along the edges of the polygons at each 
element (pixel) of the output display device. This task appropriate intersection. These calculations are some- 
becomes difficult and computationally expensive when 45 what wasteful because the values calculated are not 
the polygon to be displayed is shaded such that the themselves ultimately displayed but rather are used to 
interior of the polygon includes parameter values and is calculate other values which will be ultimately dis- 
not merely an outline of the border of the polygon. In played. This increased computational load on the com- 
this situation, the scan conversion process must fill all puter system results in slower drawing of images to a 
pixels within the polygon with appropriate parameter 50 display device. The present invention seeks to improve 
values from only the parameter values of the vertices of the method of scan converting filled polygons by avoid- 
the polygon. ing the calculation of parameter values which are not 

FIG. 1 represents a typical method of the prior art, themselves used for display, 
which is briefly discussed at pages 398-399 of Newman After a parameter value is computed for a particular 
and Sproull, Principles of Interactive Computer Graphics, 55 pixel center within a polygon, it may be necessary to 
McGraw Hill, 1979. In this technique, scan conversion modify the parameter value in order to make the ap- 
and filling begins by receiving the parameter values at pearance of the filled polygon on a display device satis- 
the vertices of the polygon, in this case a triangle as factory to a user. This is particularly true of pixel cen- 
shown in FIG. 1. For example, the vertex V0 will have ters near edges of the polygons where the original pa- 
specified with it an X and Y value (X0 and YO) and 60 rameter value of the pixel center clashes (e.g. very dif- 
parameter values such as RGB and perhaps a Z value ferent colors) with a parameter value of a pixel center 
and other parameter values. Similarly, vertex VI and which is just outside of the filled polygon. This clashing 
V2 will also have parameter values associated with the of parameter values often results in the jagged appear- 
locations of the vertices VI and V2. It will be appreci- ance of the edges of a polygon, which makes the ap- 
ated that the triangle shown in FIG. 1 is ready for scan 65 pearance of the polygon distorted on a display device of 
conversion as it has been converted to the output device the computer system. TTus appearance is largely due to 
coordinate space at subpixel resolution; thus, for exam- the fact that only a portion of the actual pixel is ob- 
ple the scan lines are shown as horizontal lines (rows) scured by the polygon and it is also due to the fact that 
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the pixel on the display device has an actual size rather 
than being an infinitely small dot. 

Prior art techniques of determining the coverage of a 
pixel and of blending or anti-aliasing the pixel within a 
filled polygon have included numerous techniques 5 
which are typically cumbersome and are often inaccu- 
rate in determining the portion of the pixel obscured by 
the polygon (i.e. in detennining coverage). It is an ob- 
ject of the present invention to provide an improved 
technique for determining the portion of the pixel ob- 10 
scured by the polygon to allow for improved blending 
and anti-aliasing of filled polygons to thereby improve 
the appearance of images displayed in the computer 
system of the invention. 

SUMMARY OF THE INVENTION 15 

The method of the present invention provides im- 
proved efficiency for scan converting a triangular poly- 
gon from parameter information at the vertices of the 
triangle in order to fill the interior of the triangle with 20 
parameter values representative of the pixel shades 
within the interior of the triangle. The method includes 
several steps, including: receiving information represen- 
tative of parameter values at each vertex of the triangle; 
selecting an edge of the triangle, which selected edge is 25 
typically a major edge of the triangle; calculating a first 
parameter value for a first pixel which is adjacent to the 
selected edge; moving to a next pixel center adjacent to 
the selected edge; calculating parameter values for the 
next pixel center adjacent to the selected edge; repeat- 30 
ing the steps of moving to a next pixel center and of 
calculating parameter values for the next pixel center 
until all parameter values for pixel centers adjacent to 
the selected edge have been calculated; and for each 
line parallel to an orthogonal axis of the display device, 35 
interpolating parameter values for pixel centers within 
the triangle which have not been calculated by interpo- 
lating from a set of calculated parameter values on each 
line. 

In a typical embodiment of the method of the present 40 
invention for scan converting filled triangular polygons, 
the selected edge is the major edge of the triangle, 
which major edge is determined by the vertices which 
specify the longest edge relative to the X axis of the 
triangle. The pixel centers which are adjacent to the 45 
major edge are selected on the basis of being as close as 
possible: to the major edge. These pixel centers which 
are adjacent to the major edge may be outside of the 
triangle, but this will not have an effect upon the ulti- 
mate result produced by the method of the present 50 
invention. 

The present invention also includes a method for 
determining the coverage of pixel centers by the poly- 
gon in order to determine the amount of blending for 
each pixel center within a polygon. The method of the 55 
invention for detennining the coverage by a polygon of 
a particular pixel center includes numerous steps, in- 
cluding detennining a value representative of the dis- 
tance from the currently sampled pixel center to each 
edge of the polygon; adding a predetermined value to 60 
each value representative of distance for the current 
pixel center; clamping each value of the distance (with 
its added predetermined value) for each currently sam- 
pled pixel center according to a function which returns 
the input of the function if the input is within a certain 65 
range and otherwise returns a value of one or zero if the 
input is out of that range; computing the coverage by 
computing the product of the clamped functions for the 
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distances; blending the parameter values if the coverage 
is within a particular range and otherwise leaving the 
parameter values unmodified for the pixel center. 

These and other aspects of the invention will be dis- 
cussed with reference to the following figures. 

DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows a triangular polygon which is over- 
layed onto an output device coordinate space at sub- 
pixel resolution for the purpose of demonstrating the 
prior art scan conversion techniques. 

FIG. 2 shows an overall flow chart of the process of 
the invention. 

FIG. 3 is a generalized architectural diagram of a 
typical computer system which may be utilized with the 
present invention. 

FIG. 4 shows an alternative computer system having 
a graphics subsystem and processor for performing scan 
conversion and other operations in connection with 
graphics processing. 

FIG. 5 shown an "up" triangular polygon as well as 
various values associate with calculations performed in 
the present invention. 

FIG. 6 shows a "down" triangular polygon which is 
being scanned converted according to the method of 
the present invention. 

FIG, 7 depicts a detailed flow chart for the scan 
conversion and filling process of the present invention. 

FIG. 8A shows a "left" triangle which is being scan 
converted according to the present invention and fur- 
ther shows the values along the major edge as well as 
interpolated values along certain rows. 

FIG. 8B shows a flowchart for the method of scan 
converting according to the present invention for scan 
conversion techniques which fill along scan lines. 

FIG. 9 shows an implementation of a hardware clamp 
of the present invention for performing' scan conver- 
sion. 

FIG. 10A conceptually shows the method of deter- 
mining the coverage of a particular pixel based upon its 
relationship with a plane. 

FIG. 10B shows the triangle of FIG. 10A on edge to 
further illustrate the concept of determining coverage 
according to the present invention. 

FIG. 11A shows a flow chart for a general method 
for determining the coverage of a pixel area by a poly- 
gon at a sample point 

FIG. 11B shows an alternative embodiment for deter- 
mining the coverage of a pixel according to the interpo- 
lation method of the present invention. 

FIG. 12 shows a triangular polygon which illustrates 
the alternative method of FIG. 11B for detennining the 
coverage of a pixel center by a polygon. 

FIG. 13 shows a quadrilateral and a pixel center 
within the quadrilateral for which the coverage of that 
pixel center is being determined according to another 
embodiment of the present invention for determining 
the coverage of a pixel by a polygon. 

FIG. 14 shows a triangle (of FIG. 10) on edge to 
illustrate the variation of the slope of the plane over the 
triangle. 

DESCRIPTION OF THE INVENTION 

The following detailed description of the present 
invention is made with reference to figures which show 
examples of the invention and which should not be 
taken as the exclusive way of performing the present 
invention. It will be appreciated by those of ordinary 
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skill in the computer art that many different implemen- gons (triangles) occurs according to the method of the 
tafaonsofthe method of the present invention are possi- present invention. After detennining the parameter 
ble and that many different computer systems may be values of each pixel center within the polygon (or alter- 
employed with the method of the present invention. natively while computing each parameter value for 
FIGS. 3 and 4 show two alternative computer systems 5 each pixel center within the polygon) the coverage by 
which may be employed with the present invention. the polygon of the particular pixel center is determined 
The computer system illustrated in FIG. 3 includes a in step 204. If this coverage indicates that the polygon 
bus 50 for communicating information, a processor 51 only partially covers the pixel then blending of the pixel 
coupled with the bus 50 for processing information, a parameter values is performed in step 205 in order to 
memory 52 such as RAM which is coupled with the bus 10 improve the appearance of the polygon and particularly 
50 for storing information and instructions for the pro- the edges of the polygon. Finally, in step 206 a display 
cessor 51, and a data storage device 54, such as a mag- file is prepared for transfer to the frame buffer which 
netic disk coupled with the bus 50 for storing informa- display file is typically a file containing the bit mapped 
tion and instructions. The computer system further image for the area covered by the polygon in the frame 
includes a video display device 53, such as a CRT raster 15 buffer. 

scan device. The scan conversion and filling of the present inven- 

Alternanve computer systems utilizing specialized tion will now be described by referring to FIGS 5 and 
graphics engines (usually specialized hardware and 6 which show an "up" triangle and "down" triangle 
software), are well known in the art. These systems respectively. As shown in FIG. 5, a triangular polygon 
typically modify the system shown in FIG. 1 by incor- 20 represented by vertices V0, VI and V2 has been con- 
porating a specialized graphics subsystem which in- verted into output device coordinate space at subpixel 
eludes a graphics processor and a dedicated frame resolution and is now ready for scan conversion accord- 
buffer, often in the form of video RAM, as well as a ing to the method of the present invention. As with the 
video display. Numerous computers from Silicon convention previously described, pixel centers are at 
Graphics, Inc. of Mountain View, Calif, exemplify 25 the intersection of scan lines and column lines. It will be 
these alternative and commercially available computer appreciated that vertex V0 has associated with it certain 
systems with specialized graphics engines. Other such parameter values, such as RGB, a Z value and perhaps 
systems are well known. FIG. 4 shows an example of a other parameter values such as an alpha value, Simi- 
computer system having a graphics subsystem. In these larly, vertices VI and V2 will also have associated 
systems, image data from the main processor 51 is com- 30 therewith similar parameter values. "P" is any one of 
municated over bus 50 and bus 60 to the graphics pro- the various parameter values (e.g., R,G,B,Z,a, etc.) P 0 
cessor 61; this information is typically in the form of represents any one of the various parameter values of 
graphic primitives such as lines, points, polygons or the vertex V0. As can be seen from the chart in FIG. 5 
character strings. The graphics processor 61 receives which shows the values along the major edge, the pa- 
nnage data from the processor 51 and uses that image 35 rameter value P 0 (which is actually any one of the vari- 
data to create an image which is rendered by using well ous parameter values at the vertex V0) is used to calcu- 
known graphic techniques such as scan conversion, late the various parameter values for each of the 8 major 
clipping, Gouraud shading and hidden surface algo- edge points which are shown in a chart of FIG. 5. The 
nthms. The resulting image is stored in video RAM 62 eight major edge points for the chart of FIG 5 are as 
(which typically includes the well known frame buffer) 40 follows: (2,4),(3,4),(4 t 4),(5,3),(6,3),(7,3),(8,2), and (9,2). 
by the graphics processor 61 addressing the video RAM The up triangle of FIG. 5 includes a major edge 
62 over the address bus 64 while supplying the video which is disposed between vertex V0 and vertex VI and 
information over bus 63. Periodically, the output of the also includes a first minor edge between vertex V0 and 
video RAM 62 is read out to a digital to analog (DAC) V2 and a second minor edge between vertex V2 and 
converter and then to a video display device or to other 45 VI. Certain values which are used in the method of the 
raster scan display devices. It will be appreciated that present invention are shown in FIG. 5; specifically, 
each scan line on the display device is parallel to an Yinit 1 and Yinit 3 are shown at points 2, 3.7 and 9, 4 1 
orthogonal axis of the display device. respectively. Yinit 1 represents the Y value at the inter- 

The computer system such as that shown in FIG. 4 is section of the major edge with the first column to the 
more appropriate for graphics intensive processing, 50 right of the vertex V0, and Yinit 3 represents the Y 
such as that required by the present invention, although value of the intersection of the first column (going from 
the present invention may operate on other computer left to right) with the second minor edge, 
systems as weU. It will be appreciated that the graphics The down triangle shown in FIG. 6 has a major edge 
processor 61 typically provides the function of scan disposed between vertices V0 and VI as well as a first 
conversion for filled polygons and thus the graphics 55 minor edge between vertices V0 and V2 and a second 
processor 61 will receive the parameter information for minor edge between vertices VI and V2. The locations 
the vertices of the polygon from the processor 51 al- of Yinit 1, Yinit2 and Yinit 3 for the triangle are shown 
though alternatively such information could also be in FIG. 6. The eight major edge points for the triangu- 
calculated by the graphics processor 61. lar polygon of FIG. 6 are shown in a chart of FIG 6 

An overview of the present invention will be de- 60 and the equations for computing the parameter values 
senbed with reference to FIG. 2. Clipped polygon in- along the major edge are also shown on FIG. 6. 
formation m the form of parameter values for each One important feature of the present invention in- 
yertex of a triangle is received by the processor which voives the observation that for each parameter (e.g. the 
is performing scan conversion (step 201). The coordi- red color value) there is a parameter planar surface 
nate values for these vertices are then converted by well 65 which may be thought of as being above the triangular 
known transformation techniques to the output device polygon because each vertex may be thought of as pro- 
coordinate space at subpixel resolution (step 202). Then viding a parameter height (e.g. the amount of red) 
in step 203, the scan conversion and filling of the poly. above a planar surface which includes the triangle in 
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that plane, which height represents the parameter value dx3 = x2-xl 

at that vertex. Since three points specify a plane, the dyl=yl-yO 

three parameter height values (for a particular parame- dy2— y2-y0 

ter, e.g. red) specify a parameter planar surface above dy3=y2-yl 

the triangle for this particular parameter. This triangu- 5 dzl=zl-zO 

lar planar surface is conceptually used to determine dz2=z2-z0 

rates of change of the parameter values through the drl=rl-rO 

triangle, which rates of change are not limited to a dr2=r2-r0 

particular edge, as in the prior art techniques* The area dgl=gl-gO 

of the triangle is actually used to calculate the rates of 10 dg2=g2-g0 

change. The psuedo-code which is presented below dbl=bl-bO 

uses this observation in order to compute the dp/dx and db2=b2-b0 

the dp/dy values (the rates of change) for each parame- dal = al-aO 

ter. da2=a2-a0 

The process of the present invention will now be 15 Note: Vertex VO at XO, YO has parameter values (rO, 

described in detail by referring to FIG. 7 along with gO, bO, zO, aO); vertex VI at XI. Yl has parameter values 

FIG. 5; this description will be for the process of the ( r l, gl, bl, zl, al); and vertex V2 at X2, Y2 has parame- 

present invention as applied to an up triangle and later ter values (r2, g2, b2, z2, a2). 

a description of the invention as applied to a down 2. compute parameter slope reciprocal 

triangle will be provided. Beginning in step 601, infor- 20 tmp=(dxl*dy2)— (dx2*dyl) 

mation for three shaded vertices of the triangle is re- (sign of tmp determines whether up or down triangle); 

ceived by the scan converting system. In step 602, the a positive tmp indicates an up triangle and a negative 

major edge is determined; often, this step is an implicit tmp indicates a down triangle (if tmp is zero, triangle 

sorting of information such that the vertices are sorted has no area and cannot be rilled) 

to provide the left most vertex of the major edge and 25 R= 1/tmp 

then the right most vertex of the major edge and finally 3. compute parameter x slopes 

the remaining vertex (labelled V2). The major edge is tmpO=dy2*R 

determined by finding the contiguous pair of vertices tmpl=dyl*R 

which has the largest difference in their values of X (if dzdx = (tmp0*dzl) — (tmpl*dz2) 

there is a tie, then arbitrarily pick one pair). Step 603 30 drdx=(tmp0*drl)— (tmpl*dr2) 

follows step 602 and it involves the determination of dgdx=(tmpO*dgl)— (tmpl*dg2) 

whether the triangle is an up or down triangle; in the dbdx=(tmpO*dbl)— (tmpl*db2) 

psuedo-code presented below, the test is performed by dadx=(tmp0*dal)— (tmpl*da2) 

detennining tie sign of the value "tmp". An up triangle 4. compute parameter y slopes 

has a positive value for trap, and the vertex (shown as 35 tmpO=dx2*R 

V2) which is not part of the major edge is above major tmpl=dxl*R 

edge. A down triangle has a negative value for tmp, and dzdy = (tmpl*dz2) — (tmp0*dzl) 

the vertex (shown as V2) which is not part of the major drdy =(tmpl*dr2) — (tmp0*drl) 

edge is below the major edge. If the triangle is neither dgdy— (tmpl*dg2)— (tmp0*dgl) 

up or down (tmp=0) then the triangle has no area and 40 dbdy=(tmpl*db2)— (tmp0*dbl) 

consequently no filing can occur. Step 603 includes a dady=(tmpl*da2)— (tmpO*dal) 

test to determine if tmp=0 and also tests if fix (x0)=fu 5. compute dy/dx along each edge 

(xl); if either of these is true, then the method is finished dydxl =dyl*(l/dxl) 

for the current triangle because no filing of that triangle dydx2 = dy2*(l/dx2) 

can occur. In step 604, the rates of change for the vari- 45 dydx3=dy3*(l/dx3) 

ous parameters are calculated for the triangle with re- Note: in this psuedo-code, derivatives are written 

spect to X and Y. Then in step 605, the parameter values without a slash, so 
for the first pixel (having coordinates at X first and Y 

first) are computed. Thus, for the example shown in df 

FIG. 5, the parameter values for the first pixel along the 50 ~3T" 
major edge are computed, which pixel is at location 2, 4; 

it should be noted that this pixel location is outside of is written as drdx; it will be understood that 
the triangle (unlike prior art methods) but that this will 

not affect the processing. This is a result of the use of dr Ar 
the parameter planar surface above the triangle which 55 ~5T* = ~Ix~ = drdK4 
will maintain a constant slope with respect to the pa- 
rameter value and hence being outside of the polygon It will also be understood that 
and using values outside of the polygon will not inter- 
fere with the computation of parameter values within dp 
the polygon. As shown in FIG. 5, the value of any one 60 "3* 
of the various parameters for the first pixel is computed 

by calculating PO+dp/dx*dxfirst-f dp/dy*dyfirst. An represents any one of drdx, dgdx dbdx, dzdx or dadx. 

example of an implementation of processing steps 603, 6. compute x of first pixel 

604 and 605 is shown below in pseudo-code which has xfirst= float (fix(x0+ALMOSTONE)) 

been annotated to further describe the invention. 65 (Note: ALMOSTONE is 0.999 . . . ) 

1. compute the delta values along the edges 7. compute y of first pixel 

dxl=xl-x0 dxfirst=xfirst-xO 

dx2 « x2-x0 yinitl =y0+ dydxl*dxfirst 
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if tmp>0 

yfirst =float (fix(yinitl + ALMOSTONE)) 
else 

yfirst= float (fu(yimtl-ALMOSTZERO)) 

Note: the value of y at the first column intersection of 

the major edge is computed in this step and is "yinitl". 
8. compute parameter values of first pixel 

dyflrst=yfirst— yO 

zfirst=zO+dxfirst*dzdx+dyfirst*dzdy 
rfirst =r0+ dxfirst*drdx + dyfirst*drdy 
gfirst = gO 4- dxfirst'dgdx + dyfirst*dgdy 
bfiret=bO+dxfirst*dbdxH-dyfirst*dbdy 
afirst = aO + dxfirst*dadx + dyfirst*dady 

Note: in the figures, a shorthand expression (e.g. 



10 



10 



15 



dp 
dy 



is shown rather than all 5 (of the last 6) equations shown 
in part 8 here. 

9. compute integer y slopes along major edge 
n-intdydxl=float (fix(dydxl)) 

10. compute two iteration slopes/parameter 
iterzl = intdydxl *dzdy -f dzdx 

iterz2 = iterzl + dzdy 
iterrl = intdy dxl*drdy + drdx 
iterr2=iterrl +drdy 
itergl=intdydxl*dgdy+dgdx 
iterg2 = itergl + dgdy 
iterbl=intdydxl*dbdy+dbdx 
iterb2=iterbl-fdbdy 
iteral — intdydxl*dady + dadx 
itera2 = iteral + dady 
Note: 



20 



25 



tterpl 



dx 



+ <«) 



Ju- 
dy 



40 



where P is any one of the parameters 2, r, g, b, or a. 

11. compute initial y of first minor edge 
yinit2 = yO + dydx2*dxfirst 

Note: the value of y at the first column intersection 
(going left to right) on the first minor edge is computed 
in this step and is "ymit2". 

12. compute initial y of second minor edge (i.e. the 
value of y at the first column intersection going from ^ 
left to right on the second minor edge is computed) 
xsecond= float (fix(x2+ALMOSTONE)) 
dxsecond = xsecond — x2 
yinit3 =y2-h dydx3*dxsecond 

The function "fix" is the integer result, rounded to 50 
minus infinity, of the input This function is used in 
several steps in the psuedo-code noted above. The 
slopes of the various parameters (e.g. RGB) with re- 
spect to X and with respect to Y are computed in steps 
3 and 4 of the psuedo-code. The computation of the 55 
parameter values for the first pixel (step 605 of FIG. 7) 
is performed in step 8 of the psuedo-code presented 
above. The value of n which used in the computation of 
parameter values along the major edge is computed in 



of "Iterpl" while the iteration value added to compute 
the pixel at 4,4 is 4 *Iterp2". The quantity "tmp" is twice 
the signed area of the triangle. Thus, for example, a 
triangle having vertices 0,4 (v0), 4,0 (vl) and 0,0 (v2) 
has a value of —16 for tmp; it can be seen that this 
triangle is a down triangle (thus a negative tmp) and has 
an area of 8. It can be seen that tmp is used to determine 
the parameter rates of change. 

Following the computation of the parameter values 
for the first pixel which is adjacent to the left most 
vertex of the major edge, the method of the invention 
proceeds to step 608 in FIG. 7 which causes processing 
to move to the next column by adding t to the value of 
xfirst. Somewhere along this next column the next pixel 
center will be found and this process is repeated for 
each column as the process "walks" along the major 
edge. This walking is performed in such a manner that 
the major edge points are as close as possible to the 
major edge in a typical embodiment and will always 
allow for filling of the triangle in an upward direction 
for an up triangle or in a downward direction for a 
down triangle. The walking along the major edge of 
FIG. 5 is shown by the arrows which jump from the 
first pixel at 2,4 (first major edge point) to the second 
pixel at 3,4 (second major edge point) to the third pixel 
at 4,4 to the fourth pixel at 5,3 and then to the fifth pixel 
at 6,3 and then to the sixth pixel at 7,3 and so on. The 
parameter values for each major edge point are com- 
30 puted as the process walks along the major edge. This 
computation is carried out in steps 617 or 618 depending 
upon whether or not the iteration value added to the 
prior parameter value is "Iterpl" or "Iterp2". Move- 
ment from major edge point to major edge point along 
the major edge and the computation of parameter val- 
ues for the major edge points continues until the current 
value of X exceeds the integer value of X for the right 
most vertex, which is XI. When this occurs, as detected 
in step 610, the walking along the major edge ends and 
the process continues by interpolating up or down the 
triangle depending on whether or not the triangle is an 
up or down triangle as determined in step 611 of FIG, 7. 
As shown in FIG. 7, an up triangle occurs when the 
variable "tmp" is greater than zero, in which case pro- 
cessing branches to node H (which has been numbered 
with the reference numeral 612); the variable tmp is 
defined in step 2 of the psuedo-code presented above. If 
the triangle is a down triangle then processing proceeds 
to node I (labeled 613). 

The logic contained within steps 609, 610, 614, 615, 
617, 618 and 619 assures that the walk along the major 
edge (i.e. the selection of major edge points) will be as 
close as possible to the major edge while remaining 
above the major edge (in the case of an up triangle) or 
below the major edge (in the case of a down triangle) 
and will assure that a major edge point can be used to 
fill all pixel centers within the triangle above the major 
edge (in the case of an up triangle) or below the major 
edge point (in the case of a down triangle). Thus, for 



35 



step 9 of the psuedo-code. The function "float" merely 60 example, the point at 5,4 in FIG. 5 would not be se- 



instructs the computer to go back to floating point arith- 
metic after deriving an integer result from the "fix'* 
function. It is noted after step 10 that for each parameter 
there are two iteration values which are computed; the 
use of these two iteration values is shown in the compu- 65 
tation of the parameter value for pixel 4,4 in FIG. 5 and 
the pixel 5,3 in FIG. 5. In particular, note that the pixel 
at 5,3 is computed by adding to P (4,4) an iteration value 



lected as a major edge point because in interpolating 
upwards from this point, the point 5,3 would be missed. 
It will appreciated that only one major edge point is 
computed for each column and that the process pro- 
ceeds left to right After the first pixel has been calcu- 
lated, initialization of certain variables occurs in step 
609. It will be appreciated that the variable YE is the 
value of Y along the major edge, which value is updated 
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as the process moves from major edge point to the next umn, as shown in FIG. 5. For example, interpolation up 
major edge point (this updating occurs in step 614). In column 3 (i.e., C3), yields: 
step 615, a test is performed to determine whether, in 

moving from the last major edge point to the next major dp 
edge point, the number of rows (scan lines) crossed by 5 p < 5 ' 4 > = Pci l = P (W + djr 

the major edge is the absolute value of n or the absolute 
value of (n + 1). It can be seen from FIG. 5 that in mov- 

ing from the first pixel at 2,4 to the second major edge _ dp 

point 3,4, the major edge does not cross a row and J>(5,5) = i>cw ~ p < 5,4 > + ~3y * 

hence the iteration value associated with "n-H 1" is used 10 

to update the value of Y and to calculate the parameter The interpolation of parameter values from major edge 
value for the pixel at 3,4 (note that n+ 1 =0 in the case points for an up triangle begins in step 625 where the 
of FIG. 5); in other words, in moving from major edge values of X and Y are reset to the first major edge point, 
point 2,4 to major edge point 3,4, step 615 selects step which is the left most major edge point. A flag (Flag v2) 

618 so that the current value of Y for the current major 15 is initialized at Not Set; this flag indicates whether pro- 
edge point remains equal to Y for the last major edge cessing of columns has gone beyond V2. 

point and the parameter values for the pixel 3,4 are Processing proceeds to step 627 in which it is deter- 
computed as shown in FIG. 5 (according to the equa- mined whether or not the first pixel (first major edge 
tion set forth in step 618). A similar result occurs in point) is within the triangle. If it is not within the trian- 
moving from major edge point at 3,4 to major edge 20 tnen processing proceeds to step 651 which incre- 
point at 4,4. After computing the new value of Y for the ments X to the next column and initializes certain height 
current major edge point and after computing the pa- (YH) and base (YB) values within the triangle, which 
rameter values for the current major edge point, the values are used to compute whether pixel centers are 
value of X is updated to move to the next column in step within the triangle. The test performed in step 627 to 

619 and the values of Y last and YE last are updated and 25 determine whether the first pixel is within the triangle 
processing loops back to step 610 which determines depends on whether X is greater than the integer value 
whether or not the walking along the major edge has of X2. If X is greater than the integer of X2 then the 
crossed beyond the right most vertex of the triangle. proper test is whether yfirst is less than or equal to 
Processing continues along the major edge through step Yinit3; in this circumstance the Flag V2 is set to its 
614 and again into step 615. It can be seen that step 615 30 " Set " state to radicate that processing from left to right 
selects step 617 when crossing from major edge point has Passed the vertex V2 (this circumstance occurs 
4,4 to major edge point 5,3. The processing continues when the lst column in the triangle is to the right of 
for all major edge points until the decision at step 610 V2 )- If x fa less than integer of X2, then the proper 
indicates that the current value of X exceeds the integer test is whether yfirst is less than or equal to Yinit 2 (in 
value of XI, which is the integer X value of the right 35 ^ V2 "« left m " Not Set " state )- The 
most vertex. The values of the points along the major appropriate test indicates a "yes" if y is less than or 
edge is described by the equations given below. e< l ual Yinit 2 X<fix(x2)) or Yinit 3 (if X>fix(x2)). It 

will be appreciated that step 627 is determining whether 

TABLE 1 the first major edge point is below the height of the up 

Equations for Values Along the Major Edge 40 triangle [defined by either the 1st minor edge (if 

X<fix(x2)) or the second minor edge (if X>fix(x2))]. 

If the first major edge point is within the triangle, 
processing proceeds to step 631 in which the parameters 
for the first major edge point are saved for filling. These 
45 parameter values may be saved in an array table or 
other means for storing the parameters. Alternatively, 
the parameters could be modified according to conven- 
tional blending techniques if the coverage determined 
by an aspect of the present invention (which is de- 
50 scribed below) indicates that some blending of the pa- 
rameter values is required. From step 631, processing 
proceeds to steps 633 and 635, which cause the process 
to move up to the next pixel center from the first pixel 
center (along the same column). Step 637 determines 
55 whether the current pixel center which was moved up 
to in step 635 is within the triangle's border; if not, then 
processing proceeds to step 651 through node K. The 
test performed in step 637 to determine whether the 
current pixel center is within the triangle depends on 
60 whether X is greater than the integer value of X2; it can 
The process of the invention continues by interpolate be seen that step 637 operates in a similar fashion as step 
ing up the triangle shown in FIG. 5 from each major 627. If the current pixel center is within the triangle's 
edge point on each column. Thus, for example, the border, then processing proceeds to step 639 where the 
major edge point at 5,3 will be used to compute parame- parameter values for this pixel center are computed, 
ter values for the pixel centers on the same column 65 The interpolation up a column according to the present 
within the triangle, which pixels centers are at 5,4 and invention merely adds the value of dp/dy to the last 
5,5. This calculation is an interpolation from the value parameter value. Then, in step 641, the parameter val- 
of the parameters at the major edge point on this col- ues for the current pixel center are saved for filling. 
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Then processing moves to the next pixel center up the 
current column by performing steps 643 and 645 and 
processing continues through step 637 until all pixel 
centers within the triangle on the current column have 
their parameter Values computed. It will appreciated 
that steps 637, 639, 641,643 and 645 are interpolating 
parameter values for pixel centers along the first col- 
umn from the first pixel (first major edge point) of the 
triangle. It can be seen from FIG. 5 that sometimes, this 
first pixel is not itself within the triangle and hence no 
interpolation up the associated column occurs; thus, for 
the triangle shown in FIG. 5 processing would proceed 
from step 625 to step 627 and then to step 651. 

The steps 651 through 703 perform the interpolation 
of parameter values from the parameter values of each 
major edge point along each column of the triangle 
(except for the first column). Step 651 calculates certain 
height related values (which indicate the height of the 
first minor edge or second minor edge at the current 
column), which values depend on whether the first 20 
minor edge or the second minor edge provides the top 
of the triangle at the current column. If the Flag V2 has 
been set, then the height is specified by the second 
minor edge, and if the Flag V2 has not been set ("Not 



14 
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15 



course be appreciated that these parameter values may 
be modified as determined by the coverage value which 
is computed as another aspect of the present invention 
(described below). After step 685, processing proceeds 
to step 687 in which the last parameter variable is up- 
dated to be equal to the current parameter values for the 
current major edge point and then the process continues 
up to the next pixel center along the current column 
(current value of X) by incrementing Y by 1 (step 689). 
In step 691, it is determined whether the current pixel 
location is within the triangle by determining whether 
the value of Y for that location is less than the height of 
the triangle. If it is, then processing proceeds to step 697 
in which the parameter values are calculated for the 
current pixel center by adding the dp/dy value to the 
parameter values for the last pixel center. It will be 
appreciated that step 697 is responsible for calculating 
parameter values of all non-major edge points up each 
column for each column after the first column, and step 
639 does this for pixels in the first column. After com- 
puting the parameter value in step 697, it may be desir- 
able to calculate the coverage by the triangle of the 
pixel center and to modify and then save those modified 
parameter values as described below or alternatively 



set" state) then the height is specified by the first minor 25 the current parameter values may be saved in step 699 



edge). Step 653 follows step 651 and determines 
whether the interpolation process should stop because 
the current column location is greater than the integer 
value of X for the right most vertex; if not, processing 
continues to step 655 in which the major edge point 30 
parameters and coordinate values are retrieved for the 
current column (X). In step 657 the process determines 
whether the current value of X has moved beyond 
vertex V2 for the first time so that the values of the 
height within the triangle need to be recomputed on the 35 
basis of the second minor edge. Thus, as shown in FIG. 
5, as the process of the invention moves from X=8 to 
X=9 the current value for YH will change as shown in 
step 657 to take into account the fact that filling up a 
column must change according to the new upper border 40 
for the triangle which is the second minor edge. Step 
657 will determine whether the current column is be- 
yond V2 for the first time by testing for the combination 
of the conditions X>fix(X2) and Flag V2=Not Set. 
Step 657 will not perform its initialization at all if the 45 
Flag V2 was set in steps 627 or 637. It is also noted that 
step 657 sets the Flag V2 if the initialization within the 
step is performed and thus this initialization will occur 
only once. In step 659, the current value for the bottom 
of the triangle is computed for the current value of X 50 
and then in step 675 it is determined whether the current 
major edge point is within the triangle. Specifically, in 
step .675 it is determined whether the current major 
edge point is below the current height for the triangle. 
Step 659 may be omitted, but the current bottom value 55 
is useful to have should it be decided to modify the 
method of the invention to allow major edge points to 
be below the major edge in the case of an up triangle; in 
this case, the current bottom value can be used to deter- 



and then processing proceeds up to the next pixel center 
on the current column through steps 701 and 703, which 
cause looping back to step 691 through node M. 

In step 691, it is determined whether the current pixel 
location is within the triangle by testing whether that 
current Y value is below the height of the triangle. If 
the current Y value is above the height of the triangle, 
then processing proceeds to step 677 which causes the 
process to move to the next column, which step also 
occurs if the major edge point for the column has a Y 
value which exceeds the height of the triangle (step 
675). If step 677 is selected, then the current X value is 
tested in step 679 to determine whether or not the cur- 
rent column is beyond the integer value of X for the 
right most vertex of the triangle, in which case process- 
ing stops; otherwise, processing proceeds to step 681 in 
which the height of the triangle is updated and then 
processing proceeds back to step 655 where the next 
major edge point is retrieved and processing continues 
for the major edge point as described above through 
steps 657 and 675. 

The interpolation process of the invention for a down 
triangle is shown in steps 751 through 825. It will be 
appreciated that these processing steps are very similar 
to the steps performed for the interpolation of an up 
triangle with certain modifications being made to the 
process in order to conform to the geometry of a down 
triangle. Specifically, the values for the base (YB) will 
have to be modified to conform with the geometry of a 
down triangle and the movement to the next pixel cen- 
ter along a particular column will occur by subtracting 
the value of one in order to move down from the major 
edge, as shown in FIG. 6. As shown in FIG, 6, X/j rr/ =2 
and Y/&tf=5; XA=dxfirst=0.9 and YA=dyfirst= -0.9; 



mine if the pixel center is below the major edge. The 60 and Y«/,i=5.6, Y to a2=2.8, and Y/„/ /3 =0.4. The coordi 



method of the invention (particularly steps 609, 614 and 
615) guarantee that the major edge points will be above 
the major edge. 

If, while on a particular column the major edge point 
is within the triangle as determined by step 675, then 65 
processing proceeds to step 685 where the parameter 
values of the major edge point for the current column 
are saved for filling of the current pixel center; it will of 



nates for Vo=(l. 1,5.9) and n= — 1. Similarly, the inter- 
polation along a particular column will merely modify 
the last parameter value by subtracting (rather than 
adding) the value of dp/dy. For instance, the parameter 
value at pixel location 5,3 shown in FIG. 6 will be com- 
puted by subtracting dp/dy from the parameter values 
at the major edge point 5,4, as shown in FIG. 6. It will 
also be appreciated that the walk along the major edge 
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and the determination of parameter values for pixel 
centers along the major edge are performed in the same 
manner as for the up triangle according to the process 
steps 608, 609, 610, 614, 615, (617 or 618), and 619. 
Thus, the major edge points of the down triangle shown 
in FIG. 6 are determined to be as follows: 
(2,5),(3,5),(4,4),(5,4), (M),(7,3),(8,3), and (9,3). The 
equations for calculating the values along the major 
edge corresponding to these points are given by Table 2 
shown below. 

TABLE 2 

Equations for Values Along the Major Edge 



10 



P(3,5) = P(2,5) + -g" + (D + 1) 



dy 



P(5,4) = 



+ 00 



dp 
dy 



+ * + 



- P(M) + -g- 



(n + 1) 



j2n_ 

dy 



p a3> = 



P (8,3) = 



P(9,3) - P(8,3) + -g- + (n + 1) 



dy 



15 



20 



25 



30 



After the computation of parameter values for major 
edge points, each major edge point for each column is 
then used to calculate parameter values for pixels cen- 
ters below the particular major edge point along a par- 
ticular column. For example, the interpolation down 
column 4 (i.e., "C4"), yields 



35 



Pom « P(3 t 3) ^5,4) - 



and 



dy 



40 



45 



It is noted that the last major edge point shown in FIG. 
6 which is at location 9,3 is actually outside of the down 
triangle; this will of course mean that no filling along 
column X=9 will occur and that this last major edge 50 
point will be discarded for purposes of filling the trian- 
gle. 

A feature of the present invention allows for the 
computation of parameter values outside of the triangle 
to be used for computing parameter values within the 55 
triangle. FIG. 5 shows an example of this where the first 
pixel center is outside of the triangle and its parameter 
values are used to compute all other parameter values. 
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order to prevent clamping by floating point systems in a 
computer system. 

It will be appreciated that certain types of specialized 
graphics hardware may be used with the present inven- 
tion. For example, parallel processing of each column 
separately by separate processors (after the calculation 
of parameter values for a major edge point on a column) 
may be used in order to more rapidly calculate the 
interpolated parameter values along each column. 
However, the method of the present invention may be 
performed with general purpose computer systems and 
does not require such specialized graphics processors or 
parallel processing. It will be appreciated that numerous 
modification to the method of the present invention are 
possible. For example, parameter values up or down a 
column (depending on whether the triangle is an up or 
down triangle) may be computed immediately after 
computing the parameter values for the major edge 
point on the column (rather than waiting for the calcu- 
lation of parameter values for all major edge points). 
The implementation of this approach is readily within 
the capabilities of those with ordinary skill in the art 
upon reference to the method disclosed herein. 

The method of the invention should not allow clamp- 
ing of parameter values during interpolation, as de- 
scribed above, and specialized hardware may be used 
with the method of the present invention to assure that 
after obtaining all interpolated parameter values, these 
values are checked to detennined whether they have 
gone beyond maxima or minima allowable values. An 
example of an implementation of this hardware for 
detenriining whether allowable maxima or minima have 
been exceeded is shown in FIG. 9. The register 90 is a 
multibit register which receives the output of each pa- 
rameter value at pixel centers within the triangle from 
the triangle filling process of the present invention. The 
interpolation process uses an extra bit beyond the ex- 
pected most significant bit (for each parameter value) 
and this extra bit is provided to the register 90 along 
with the rest of the parameter value. It will be appreci- 
ated that several such registers (and the other associated 
logic shown in FIG. 9) may be implemented for parallel 
receipt of the various parameter values for a pixel cen- 
ter. The parameter value is stored as a multibit value in 
register 90, where the most significant bit of the parame- 
ter value is stored in the bit register 91 which is part of 
the register 90. Register 90 also includes an overflow bit 
register 92 which will indicate whether or not the pa- 
rameter value has overflowed beyond the allowable 
maxima or minima. The value of the overflow indicator 
bit in register 92 provides a signal over the select line 94 
of the multiplexor 95, which signal selects between the 
two inputs from buses 93 or 98 to the multiplexor 95. 
Specifically, if the value of the overflow indicator bit in 
register 92 is a zero, then the parameter value, including 
the most significant bit (MSB) of the parameter value in 
register 91, is provided over bus 93 and out of multi- 
plexor 95 via bus 101. In other words, if the overflow 



As explained above, because of the "planar parameter 

surface" located conceptually above the triangle, this 60 indicator bit in register 92 is a zero, then the parameter 

computation will still remain accurate in going from value is selected by the multiplexor 95 for output. If on 

parameter values outside of the triangle to parameter the other hand, the overflow indicator bit is a one, then 

values within the triangle. It should be noted that dur- the value selected by multiplexor 97 is provided at the 

ing this interpolation process, parameter values should output of multiplexor 95. The output of multiplexor 97 

not be clamped (i.e. they should be allowed to fluctuate 65 is selected on the basis of the status of the select line 96 

over as wide a range as possible with the computer which is coupled to receive the value of the most signifi- 

system being used). Preferably, fixed point arithmetic cant bit of the parameter value as stored in location 91. 

may be used during the process of the interpolation in If the most significant bit is a one, then the input 99 to 
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multiplexor 97 is selected for output by MUX 97. Con- putations. Thus, for example, the equation for "tmp" for 

versely, if the most significant bit in location 91 is a a row triangle is: 
zero, then the input 100 is selected as the output of 
multiplexor 97. As shown in FIG. 9, the input 99 is all 

zeros while the input 100 is all ones. It will be appreci- 5 tmp = (^0(^2) - «f2)(^i) 

ated that if the most significant bit is a one when an = 0>1 - X))(x2 - *0) - o* - j€)(*i - x0) 

overflow condition occurs then the proper value of the A1 , , , , 

parameter is most likely all zeros, which are selected by M ™ shown below are examples of equations for a 

the multiplexor 97 and then by the multiplexor 95 for row tnan * wl ?ch we 5 e den t ved ^ substitution of 

output from multiplexor 95. If an overflow condition » *** J £ *) *om the corresponding column 

occurs when the most significant bit is a zero, then the tm^O^dxZ^R 00 

most likely parameter value will be all ones (input 100 tmDl==dxl*R 

to multiplexor 97) which is provided as the output of dr d y =(tmp0*drl)-(tmpl*dr2) 

multiplexor 95. drdx=(tm P l*dr2)-(tmpO*drl) 

An example of the scan converting and filling process dxdyl=(dxl)*(I/dyl); dxdy2==dy2(l/dy2) ; 

of the present invention in the context of triangles dxdy3=dx3(l/dy3V 

which are filled along scan lines ("rows**) rather than yfirst= float (fix(y0-t-O99 )) 

columns will now be described with reference to FIGS. dyfirst=yfirst-yO 

8A and 8B. Triangles filled along rows are referred to as 2Q xinitl=x0+dxdyl*dyfirst 

"row triangles" and triangles filled along columns are if i e ft triangle (tmp<0): xfirst afloat (fix(xinitl 

referred to as "column triangles." FIG. 8A shows a — ALMOSTZERO)) 

'left" triangle having a major edge between vertices V0 if right triangle (tmp>0): xfirst=float (fix(xinitl+0.99 . 

and VI and a second minor edge between vertices V2 . . )) 

and VI and a first minor edge between vertices V2 and 25 dxfirst = xflrst — xO 

V0. It can be seen that X/? /5/ =9 and Y/?„,= 2. According n=float(fix(dxdyl)) 

to FIG. 8A: dxfirst=XA= -0.8, dyflrst=YA=0.9, xinit2 = xO -f dxdy2(d yfirst) 

Xinit3==x2+dxdy3[float(fix(y2+0,99 . . . ))-y2] 
It will also be appreciated that the equations used in 

„ „ Fixf dx\ \ _ lf 30 the various steps of FIG. 7 will be similarly modified 

\ dyX ) ' (substituting x for y and y for x) in order to account for 

the geometry of a row triangle. Since 

and n=Fix(dx dyl). The first pixel location, represent- 
ing the first major edge point is at 9,2 and there are four dp 
other major edge points which are determined in se- 35 
quence from the first major edge point. These four 

major edge points are: (9,3),(9,4), (8,5), and (8,<5). The subtracted from parameter values in interpolating left 

major edge is determined by finding the pair of vertices dcm S a row in a left triangle, a left triangle is analogous 

with the largest difference in their values of Y (in the to a down trian g! e (where 
case of a tie, pick one pair). The triangle of FIG. 8 A is 40 

a left triangle because the vertex (shown as V2) which ji 
is not part of the pair of vertices specifying the major 

edge is to the left of the major edge. The values along . „ . . . , 

the major edge for the major edge points are shown in 18 2^?^ ******** down a column . 

FIG. 8A and are calculated in essentially the same man- 45 The mterpolatxon along each scan line (which is an 

ner as the parameter values along the major edge points ortho S™* ™* <* ^ display device) occurs after the 

for the do™ triangle discussed previously. Theequa- Parameter values for the major edge point on that scan 

tions are given in TWe 3 below *? e ^ ve b ? n ^P^ed In the case of the left mangle, 

^ the interpolation proceeds by subtracting dp/dx from 

TABLE 3 50 the last parameter value. For instance, the calculation of 

Equations for Values Along the Major Edge the parameter values are given as 



dy 



p<«) 


= P0WUTO-l--i£-X 4 + -ijE-Y A 


P <M) 


= Pp.2) + -g- + <n-H)-g- 


P <W) 


= P(W) + -^- + (° + 0-^- 


*<M> 






= *<8.S) + + 



55 



Prs-\ = P<8,6) - ; Prl\ = P(s,s) - -fe- ; 



and 

60 PR^^PRV-gr 

FIG. 8B generally illustrates the method of the inven- 
tion for triangles having parameter values which are 
interpolated along scan lines. The parameter values for 
In practice, the computations for row triangles are the 65 the triangle's vertices are received in step 500 and then 
mirror image of computations for column triangles be* the major edge of the triangle is determined in step 501, 
cause "x w is substituted for all *y s" in column computa- as well as the determination of whether the triangle is a 
tions and **y M is substituted for all "xV* in column com- left triangle or a right triangle. A right triangle has a 
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positive value for tap, and the vertex which is not part (pixel center) 1003 is also two (therefore, pixel center 
of the major edge is to the right of the major edge. A 1003 is two pixels into the interior of triangle 1001 ). 
left triangle has a negative value for tmp, and the vertex FIG. 10B shows 0 another view of the triangle 1000 and 
(shown as V2) which is not part of the major edge is to the plane 1002. This view is taken from on the edge of 
the left of the major edge. In step 502, the parameter 5 triangle 1000 looking down edge 1007. From this view, 
values for the first pixel center adjacent to the major edges 1001 and 1005 appear as points and the plane 1002 
edge are computed. In the case of the left triangle appears as a line. It can be seen from this view that the 
shown in FIG. 8A, the value at location 9,2 is calculated height between the two planes and the distance of the 
according to the equation shown in FIG. 8A. Then, in pixel from the edge will be equivalent because the slope 
step 503, the process moves up to the next scan line and 10 is equal to one. The height will have a negative value if 
the process picks the next pixel center along the major measured from a point outside of the triangle. This 
edge; this continues until each scan line intersecting the height is an approximation of the coverage by the poly- 
major edge is provided with a major edge point along gon of the particular pixel center and it is a value repre- 
the scan line. At each new major edge point, the param- sentative of the distance of the pixel center from the 
eter values are calculated in step 504 according to the 15 edge which intersects the plane. It will be appreciated 
method of the present invention, utilizing the equations that a higher height indicates that it will be less likely to 
shown in FIG: 8A. Then processing proceeds to step filter the parameter value of the pixel center. 
505 where parameter values are interpolated along each One method of detennining this height is to compute 
scan line containing a major edge point. The interpola- and sample at each pixel center which is to be filled two 
tion of parameter values along each scan line containing 20 plane equations, one representing the plane 1002 and the 
a major edge point proceeds in a normal fashion, includ- other representing the plane 1010, which is the plane of 
ing a test to determine whether the pixel center is within the triangle 1000. These two planes will have a slope 
the triangle. It will be appreciated that scan line inter- with respect to each other of one. For each pixel center 
sections with the edges will be utilized to determine the within the polygon, the distance between the two 
m a xim u m left and right values for the interpolation 25 planes is computed and this provides a value for each 
process and to test if pixel centers are within the trian- pixel center which is proportional to (and a function of) 
gle. The testing step 507 determines whether the last the distance from the pixel center to the edge which 
row has had its pixel centers interpolated and if not intersects the plane 1002. This step is repeated for each 
processing loops back to step 505 until all scan lines edge so that each pixel center has a measure of distance 
having pixel centers within the triangle are interpolated. 30 relative to each edge of the polygon. These measures of 

After completing all interpolations of pixel centers distance are then used according to the process of the 
within the triangular polygon, the values for the param- present invention to provide an estimate of the coverage 
eters of all pixel centers within the polygon may then be by the polygon of the particular pixel. Other methods 
assembled in a display file for transfer to the frame for detennining the distance are discussed below, 
buffer where conventional operations may be per- 35 A general method of the present invention for deter- 
formed, such as hidden surface algorithms utilising the mining the coverage by a polygon of a pixel will be 
Z values and other well known graphics techniques. described while referring to FIG. 11A. It will be appre- 
Alternatively, the parameter values may be filtered dated that the method of the present invention for de- 
according to conventional anti-aliasing techniques in tenninicg coverage is applicable to triangular polygons 
order to improve the appearance of the polygons, par- 40 as well as nontriangular polygons, such as the quadrilat- 
ticularly the edges of the polygons. Thus, as shown in eral polygon shown in FIG. 13. The method begins in 
FIG. 2, after completing the scan conversion and filling step 1101 which may be either after filling pixel centers 
process of the present invention, in step 203, the param- within the polygon or during the process of filling pixel 
eter values may be sent to the frame buffer for display centers within the polygon. In step 1102, the process 
or alternatively filtering of the parameters values may 45 determines a value representative of (and proportional 
be performed based upon the determination of the cov- to) the distance from the currently sampled pixel center 
erage of each pixel center by the particular polygon. to each edge of the polygon. The distance dl represents 
These steps are illustrated in steps 204 and 205. An a function of the distance to edge 1, the distance d2 
aspect of the present invention includes an improved represents a function of the distance to edge 2, and the 
method for detennining the coverage by a polygon of a 50 distance dn represents a function of the distance to edge 
particular pixel and this technique may be included into n. These distance values may be proportional to (but not 
the interpolation method described above so that cover- equal to) the actual distance, depending on the slope of 
age is determined for each pixel center within the poly- the plane which intersects the edge. If the slope is one, 
gon during the filling process. then the height will equal the actual distance to the edge 

The method of determining the coverage of a pixel by 55 (as shown in FIG. 10B). However, if the height value 

a polygon will now be described by referring to FIGS. is computed according to a plane (which intersects the 

10A and 10B. Conceptually, imagine the plane of a edge to which the distance value is being computed) 

triangle being intersected at one of its edges by a plane which has a slope of other than 1, then the distance 

having a slope of one relative to the plane of the trian- value represented by the computed height will be pro- 

gle. This is shown in figure 10A wherein the triangle 60 portional to (but not equal to) the actual distance to the 

1000 is intersected at its edge 1001 by the plane 1002. edge. For example, in FIG. 14 a plane 1020 (which 

Plane 1002 includes the edge 1001 and has an angle with intersects edge 1001) has a slope of two relative to the 

respect to the plane of the triangle of 45 degrees (equiv- plane 1010 which is the plane of the triangle 1000. In 

alent to a slope of one). Point 1003 is within the plane of this case, the distance value represented by the com- 

the triangle and is a pixel center within the triangle. 65 puted height will be twice the actual distance (from the 

Point 1004 is a point in the plane 1002. It is clear that the point in the triangle to the edge), and the edge will 

"height" from the point 1003 to the point 1004 is two appear (after blending) sharper. In the preferred era- 

and that the distance from the edge 1001 to the point bodiment, the slope of the plane (which intersects the 
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edge of the polygon) is 1. If a slope less than 1.0 is used, 
then more blurry filtering (blending) is achieved. In step 
1103, a predetermined value is added to each value 
representing distance for the currently sampled pixel 
center. In the currently preferred embodiment of the 5 
present invention,the predetennined value is 0.0, al- 
though other values (e.g. 0.5, etc.) may be chosen. The 
predetermined value which is added to the value repre- 
senting distance in step 1103 is added to handle situa- 
tions where pixel centers are close to the edge. Pixel 10 
centers close to an edge (e,g. pixel center at 5,3 of FIG. 
5) will have height values of less than 1 and thus will 
tend to produce clamped functions of less than one 
(resulting in more blending at these pixel centers); add- 
ing the predetermined value tends to increase the cover- 15 
age of these pixel centers. A predetennined value which 
is closer to zero will tend to shrink the polygon. In step 
1105, each value representative of distance (as modified 
by step 1103) is then clamped for the currently sampled 
pixel center according to a function which clamps the 20 
output of the function to a predetermined value (e.g. 1) 
if the input to the function is outside of a predetermined 
range. In the preferred embodiment, the value represen- 
tative of the distance is clamped to the value of one if 
the input to the function is greater than one. Also in the 25 
preferred embodiment, the value of the clamp function 
is zero if the value of the input to the function is less 
than zero. In the preferred embodiment, the value of the 
clamp function is equal to the input to the function if the 
input is within a range of zero to one, inclusive. 30 

Following the clamping step in 1105, the value of the 
coverage is computed in step 1106 by multiplying all 
output values of each clamped distance function. Then 
in step 1107, the parameter values (e.g. RGB) are 
blended according to conventional prior art techniques 35 
if the coverage is within a predetermined range and 
otherwise the polygon is filled with the currently calcu- 
lated parameter values which remain unblended. In the 
currently preferred embodiment, the predetermined 
range for blending is: 0<coverage<l. In other words, 40 
blending will occur if the value of coverage is between 
0 and 1 (not including the extremes of that range). If the 
value for coverage is equal to 0 then no filling of the 
pixel will occur and if the value of the coverage is equal 
to 1, then the polygon is filled at that pixel location with 45 
the currently calculated parameter value which remains 
unblended. Blending is performed using well known 
pixel blending techniques which use, as an input, the 
value of the coverage of a pixeL For example, a simple 
blending technique for a gray scale display device may 50 
involve multiplying the gray scale parameter for a pixel 
by the coverage of the pixel. Thus, if the coverage =0.5, 
then the blended pixel value is one-half of the original, 
unblended gray scale value. 

A specific method of the present invention of deter- 55 
naming pixel coverage will now be described by refer- 
ring to FIGS. 11B and 12. This method performs the 
computation of coverage and blending during the inter- 
polation process of the present invention and utilizes the 
same interpolation techniques used for scan converting 60 
and filling triangular polygons. The three distances of a 
point to the three edges of the triangle are treated as 
three separate parameters which are interpolated ac- 
cording to the method shown in FIG. 7 above. FIG. 12 
shows an up triangle having vertices V0, VI and V2. 65 
Vertex V0 has the coordinate location of X0, Y0, vertex 
VI has the coordinates XI, Yl, and vertex V2 has coor- 
dinates X2, Y2. It is defined that Xf mt is 2, Y/fnr is 4, 
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XA=dxfirst=0.3, YA=dyfirst=0.2, and n= — 1. In 
addition, it is given that A at V0=0; B at V0=Bo; C at 
V0=0; A at V1=0; A at V1=0; B at VI =0; C at 
Vl=Co; A at V2=Ao; B at V2=0; and .C at V2=0. 
This method of the present invention is based upon the 
observation that the plane which intersects an edge of 
the triangle, which plane has a slope of one relative to 
the plane of the triangle, will have a height (above the 
triangle at the opposing vertex) which is equal to the 
distance from the opposing vertex to the edge. Thus, 
each vertex has a value representative of the height of 
the three planes above that vertex. Thus, for example, 
the plane having a slope one and intersecting edge B 
will have a height above vertex VO of B0 (while verti- 
ces VI and V2 have a height equal to zero relative to 
this plane). It can be seen that this plane is associated 
with the parameter value B which represents the func- 
tion of the distance from a point to tie edge B; notice 
that vertices VI and V2 have a distance of zero. As 
shown in FIG. 12, the parameter A represents the dis- 
tance from a point to edge A, B represents the distance 
from a point to edge B, and the parameter C represents 
the distance from a point to edge C. Bo is defined as the 
distance from vertex V<j to edge B. The value of BO may 
be derived in various ways, including the method de- 
scribed previously of computing the two plane equa- 
tions and then determining the distance between vertex 
VO in the plane of the triangle and the corresponding 
point in the plane which intersects edge B. An alterna- 
tive method for computing B0 is the method described 
in conjunction with FIG. 13, which computes the inter- 
section on edge B of the perpendicular line containing 
the vertex VO, which line is perpendicular with respect 
to the edge B. The three parameter values A, B and C 
are shown at the three vertices in FIG. 12; AO is the 
distance from V2 to edge A, and CO is the distance from 
VI to edge C. 

FIG. 11B shows the method of determining pixel 
coverage during the interpolation process of the present 
invention. As can be seen from FIG. 12, the interpola- 
tion of pixel coverage values is performed in the manner 
shown in FIG. 7 except that the parameter slopes for A, 
B and C with respect to X and Y are computed as 
shown in FIG. 12 (rather than as shown in conjunction 
with steps 604 and 605 of FIG. 7). Moreover, in step 
1152 of FIG. 11B, the computation of Afirst, Bfirst and 
Cfirst (for the first major edge point) is performed as 
shown in FIG. 12 (rather than as shown in step 605 and 
in FIG. 5). The process of FIG. 11B is executed each 
time parameter values are computed in the method of 
FIG. 7, including steps 605, 617, 618, 639, 697, 763 and 
819. Table 4 lists the equations describing the A,B, and 
C parameters values for the major edge points. 

TABLE 4: Equations for Values Along the Major Edge 

A(i t i)-A(4 t 4)-\-dAdx+{n)dAdy 
^3)=^(5,3)+<W^H-(«+ \)dAdy 
^.3)=^(6.3)+ dAdx+ („+ \)dAdy 
A$2)=AQ2)+dAdx+{n)dAdy 
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Bftrst = 5(2,4) — -^0 + (dBdx)X& 4- (dBdy) Y& 
B$ t <i)=Bw\*dBdx+(n+ \)dBdy 
5(4,4)=5(3,4)+^^+(n+ 1)d£aTK 

5(7,3)= ^(6,3) + \)dBdy 

Cfint= C^^{dCdx)X^{dCdy)Y L 
c (M)=C(2,4)+^C<«x+(n+ lyafy 
%4)=C[3,4)+rfC^+(n+ \)dCdy 
CwaCptf+dCdx+WdCdy 
C(6,3)=C( 5t 3)+rfCd!r+(n+ l)dCrfy 
07,3)- C m +dCdx+("+ WCdy 
C(8 t 2)=C(7,3)+^C^+(nyCrfjr 
^9 f 2)=C(8,2)+</Cax+(rt+ l)dCDtK 
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The parameters for the pixel centers are determined 
through interpolation. For example, interpolating up 
column X=5 yields: A(5,4)=A(5,3)+dAdy; A^s. 
)« A( 5 , 4 )+dAdy; B (5 ,4)«B( 5> 3)H-dBdy; 35 

B (5,5)=B( 5 ,4)+dBdy; q 5 ,4)=C ( 5 ) 3)+dCdy; and 
C( 5t 5)«=C(5,4)+dCdy. 

After computing the A, B and C parameter values for 
each pixel center, the following distance parameter 
slope equations are utilized. 40 



dAdx = 



dBdx = 
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dCdx-- 
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where 



dx\ = XX - X0 dy\ = Y\ - 3D 
dxl = XI - X0 dyl = n - Y0 
dx2 - X2 - X\ dyl = Yl - Y\ 
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proceeds to step 1153 in which the three parameters A, 
B and C may be modified according to the slope of the 
plane intersecting each edge, as shown in step 1153. In 
the preferred embodiment, the slope is equal to one and 
thus there is no modification of these parameters. In step 60 
1154, a predetermined value is added to each of the 
three parameters for each pixel center. In the preferred 
embodiment, this predetermined value is 0.0. Then in 
step 1155, the clamp function of each of the three pa- 
rameters A, B and C is computed for each pixel center 55 
according to the clamp function specified in step 1155. 
In particular, the clamp function provides the value of 
the input if the input is within a predetermined range, 



which is typically zero to one, inclusive. The clamp 
function returns the value of zero if the input is less than 
zero, and the clamp function returns the value of one if 
the input is greater than one. 

In step 1156, the coverage at each pixel center is 
computed by multiplying the clamp function of each of 
the three parameters A, B and C for the particular pixel 
center. Finally, parameter values are blended if the 
coverage falls within a predetermined range as shown 
in step 1157. 

The method of determining pixel coverage by a poly- 
gon, including an N sided polygon, where N exceeds 
three, may be described generally with reference to 
FIG. 13. FIG. 13 shows a quadrilateral polygon having 
vertices 1251, 1252, 1253 and 1254 which form edges 
1255, 1256, 1257 and 1258. This polygon 1250 is shown 
in an output device coordinate space having a scan line 
at YC, which scan line has been labeled 1261, The inter- 
section of the current X column with the scan line 1261 
finds the pixel center XC, YC for location 1260. A 
method will now be described for the general computa- 
tion of the distance from a pixel center within a polygon 
to any edge of the polygon. This distance may then be 
used as part of the process shown in FIG. 11A to deter- 
mine the coverage by the polygon of the particular 
pixel. This method attempts to find the line through the 
current pixel center to the particular edge, which line is 
perpendicular with respect to that edge. It is known that 
such a line will have a slope which is the negative recip- 
rocal of the slope for the particular edge. If the slope of 
the edge is zero or infinity (division by zero) then the 
calculation is simplified. If the slope of the edge is zero 
then the line of the edge is horizontal and therefore the 
distance is equal to the absolute value of the quantity (Y 
value at one of the vertices of the edge minus the Y 
value at the current pixel center). If the edge has an 
infinite slope then the edge is vertical and the distance 
to the edge from the current pixel center is simply the 
absolute value of the quantity (X of one of the vertices 
minus XC). If neither of these situations is the case, then 
the value for B (from the equation Y=MX+B) for the 
line which perpendicularly intersects the edge is com- 
puted, as follows: 

B c =McK c -Yc 

It will be appreciated that XC and YC are known and 
that MC is the negative reciprocal of the slope of the 
edge to which the distance is being calculated. After 
determining Be, then two equations are simultaneously 
solved for the value of the intersection point (XI, YI) on 
the edge with the perpendicular line. Specifically, the 
simultaneous solution of Y=McX+BC and Y=M^ 
X-fB£ is calculated (where ME is the slope of the 
edge). It will be appreciated that 



where Xn and Yn are the X and Y coordinates at one 
vertex of the edge and Xvo and Yvo are the coordinates 
at the other vertex of the edge. The simultaneous solu- 
tion gives the intersection point (XI, YI) of the perpen- 
dicular line with the edge which then permits the com- 
putation of the distance according to the standard Eu- 
clidean distance formula— distance = ((XI — XC) 2 - 
+(Y*— YC) 2 )i. This computation is performed for each 
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perpendicular intersection with each edge from the 
point XC, YC. After all distances from that point to 
each edge are computed, then the method of the inven- 
tion may be utilized to determine the coverage of the 
pixel at XC, YC by the polygon by performing steps 
1103, 1105, and 1106 to determine the value of cover- 
age. Then, in step 1107, conventional blending tech- 
niques may be utilized to filter the parameter values at 
the current pixel location. 

It will be appreciated that other well known tech- 
niques for computing the distance of a pixel center to 
the edge of a polygon may be utilized in accordance 
with the present invention for determining the coverage 
of a pixel by a polygon. It will also be appreciated that 
other equivalent processing steps in the method for scan 
converting and filling a triangular polygon maybe sub- 
stituted for the specific steps presented in FIG. 7. It 
should also be understood that the figures are for pur- 
poses of illustration only and should not be taken as 
limitations upon the invention. It is clear that the meth- 
ods and apparatus of the present invention may be mod- 
ified and changed in accordance with uses desired from 
the present invention. It is contemplated that many 
changes and modifications may be made by one of ordi- 
nary skill in the art without departing from the spirit 
and the scope of the invention. 

I claim: 

1. In a computer controlled display system with a 
display device having an array of pixels, a method for 
rendering a polygon for display said method comprising 
the steps of: 

a) deterniining a coverage of a pixel within said poly- 
gon by said polygon by: 

i) determining a first value representative of a dis- 
tance from a center of said pixel to an edge of 35 
said polygon; 

ii) adding a first predetermined value to said first 
value representative of the distance for said pixel 
to produce a second value.; 

iii) clamping said second value according to a tunc- 40 
tion which provides a constant second predeter- 
mined value if said second value is less than a 
first predetermined range and which provides 
the second value if the second value is within 
said first predetermined range and which pro- 45 
vides a constant third predetermined value if said 
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second value is greater than said first predeter- 
mined range; 
iv) repeating steps i— iii for each edge of said poly- 
gon; 

5 v) computing the coverage of the pixel by calculat- 
ing the product of all clamped values corre- 
sponding to edges of said polygon; 
b) determining coverages of each pixel of said poly- 
gon; 

10 c) storing a plurality pixel values in memory, wherein 
said pixel values are a function of said coverage; 

d) generating electrical signals to actuate pixels of 
said polygon according to said pixel values; 

e) displaying said polygon on a display device, 
15 2. A method as in claim 1 wherein said first predeter- 
mined value is 0.0 and wherein said second predeter- 
mined value is 1.0. 

3. A method as in claim 2 wherein said first predeter- 
mined range is from zero to 1, inclusive. 
20 4. A method as in claim 3 further comprising the step 
of blending parameter values for a pixel if the value of 
said coverage is within a second predetermined range. 

5. A method as in claim 4 wherein said second prede- 
termined range is from zero to 1, not including the 

25 extremes of said second range. 

6. A method as in claim 1 wherein said step of deter- 
mining a first value representative of a distance from a 
center of said pixel to an edge of said polygon is per- 
formed during the scan conversion and filling process 

30 for said polygon. 

7. A method as in claim 1 wherein the step of deter- 
niining a first value representative of a distance from a 
center of said pixel to an edge of said polygon is per- 
formed by determining a height projecting away from 
said polygon at each pixel center of a point directly 
above said each pixel center on a plane intersecting each 
edge of said polygon. 

8. A method as in claim 7, wherein said plane has a 
slope of 1 relative to the plane of said polygon. 

9. A method as in claim 7, wherein said polygon is a 
triangle. 

10. A method as in claim 1 wherein said first value 
representative of a distance from a center of said pixel to 
an edge of said polygon is proportional to said distance. 

***** 
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